Estimation of the drift of fractional Brownian motion 
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Abstract 

We consider the problem of efficient estimation for the drift of fractional 
Brownian motion B H := (-B t ff ) te ^ T j with hurst parameter H less than i. We 
also construct superefficient James-Stein type estimators which dominate, under 
the usual quadratic risk, the natural maximum likelihood estimator. 
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1 Introduction 

Fix H e (0, 1) and T > 0. Let B H = Ub? '\ flf 4 ); t G [0, T]\ be a d-dimensional 

fractional Brownian motion (fBm) defined on the probability space (£l,J-,P). Lhat 
is, B H is a zero mean Gaussian vector whose components are independent one- 
dimensional fractional Brownian motions with Hurst parameter H £ (0,1), i.e., for 
every i = l,...,d B H ' 1 is a Gaussian process and covariance function given by 

E(B^B^) = I (s 2H + t 2H -\t- s\ 2H ) , s,te [0,T\. 

For each i = 1, ... ,d, {^t) t e[0T} d en °t es the filtration generated by (^Bf' l ^j ^ . 
The fBm was first introduced by [S] and studied by [B]. Notice that if H = i, the 
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process Bz is & standard Brownian motion. However, for H ^ i, the fBm is neither 
a Markov process, nor a semi-martingale. 

Let M be a subspace of the Cameron-Martin space defined by 

M = j</? : [0, T] -> = jf #ds with ^ G L 2 ([0, T]) 

and ^€i^ + ' (L 2 ([0,T])),i=l,...,d\. 

Let 6> = {($,... ,0f);t G [0,T]} be a function belonging to M. Then, Applying 
Girsanov theorem (see Theorem 2 in [9] ) , there exist a probability measure absolutely 
continuous with respect to P under which the process B H defined by 



B* I = B?-e t , te[o,T) (i.i) 

is a fBm with Hurst parameter H and mean zero. In this case, we say that, under 
the probability Pg, the process B H is a fBm with drift 9. 

We consider in this paper the problem of estimating the drift 9 of B H under 
the probability Pg, with hurst parameter H < 1/2. We wish to estimate 9 under the 
usual quadratic risk, that is defined for any estimator 5 of 9 by 



K(9,5)=E d 



T 

\\5 t -e t \\ 2 dt 







where Eg is the expectation with respect to a probability Pg. 

Let X = (X 1 , ... , X d ) be a normal vector with mean 9 = (0 1 , . . . , 9 d ) G R d 
and identity covariance matrix a 2 Id- The usual maximum likelihood estimator of 9 
is X itself. Moreover, it is efficient in the sense that the Cramer-Rao bound over all 
unbiased estimators is attained by X. That is 

a 2 d = E\\\X-e\\ 2 d ] =ME[u-e\^[, 

5fc<-> 

where S is the class of unbiased estimators of 9 and denotes the Euclidean norm 
on R d . 

constructed biased superefficient estimators of 9 of the form 
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for a sufficiently small and b sufficiently large when d > 3. [4 J sharpened later this 
result and presented an explicit class of biased superefficient estimators of the form 



i 1 ~ pJp) X ' f ° r ° < a < 2 ^ ~ ^' 



Recently, an infinite-dimensional extension of this result has been given by 
[TO] , The authors constructed unbiased estimators of the drift (Ot)t£[oT] °f a con " 
tinuous Gaussian martingale (Xt) t£ ^ T ^ with quadratic variation a^dt, where a £ 
L 2 ([0, T], dt) is an a.e. non- vanishing function. More precisely, they proved that 
9 = (X t ) te [ 0T j is an efficient estimator of (0t)te[0T]- On the other hand, using Malli- 
avin calculus, they constructed superefficient estimators for the drift of a Gaussian 
processe of the form: 



X t := I K(t,s)dW s , t€ [0,T], 
J o 



n(9,B H )=Ee 



-d. 



where iWt)t^[o,T] 1S a standard Brownian motion and K(., .) is a deterministic kernel. 
These estimators are biased and of the form Xt + DflogF, where F is a positive 
superharmonic random variable and D is the Malliavin derivative. 

In Section 3, we prove, using technic based on the fractional calculus and 
Girsanov theorem, that 6 = B is an efficient estimator of 9 under the probability 
Pq with risk 

r r T 1 rp2H+l 

/ \\B?-6 t fdt =±= . 

Jo " * " J 2H + 1 

Moreover, we will establish that 9 = B H is a maximum likelihood estimator of 9. 

In Section 4, we construct a class of biased estimators of James-Stein type of 
the form 

5(B% = (l - at™ ( ! ^{?)) Bf, t E [0,T]. 

We give sufficient conditions on the function r and on the constant a in order that 
5(B H ) dominates B H under the usual quadratic risk i.e. 

TZ(9,6 (B H )) <TZ(9,B H ) for all 9 £ M. (1.2) 



2 Preliminaries 

This section contains the elements from fractional calculus that we will need in the 
paper. 
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The fractional Brownian motion B has the following stochastic integral representa- 
tion (see for instance, [I], [8]) 



B?'* = [ K H {t,s)dW l s , i = l,...,d; t € [0,T\ (2.3) 
Jo 

where W = (W 1 , W d ) denotes the d-dimensional Brownian motion and the kernel 
Kn(t, s) is equal to 

c H (t - s) H - l i + c H { 1 - - H) f\u - s) H ~'i (l - ( du if H < ~ 

c H {H - 1) J\u - s) H ~i du i£H>± 



if s < t and Kn(t, s) = if s > t. Here c# is the normalizing constant 



/ 2flT(§ - H) 
T(H+±)T(2-2H) 



where T is the Euler function. 

We recall some basic definitions and results on classical fractional calculus which we 
will need. General information about fractional calculus can be found in 
The left fractional Riemann-Liouville integral of / G L 1 ((a, 6)) of order a > on 
(a, b) is given at almost all x G (a, 6) by 

1 f' x 

1 1 ' \a-l 



^ +f{x) = f{a)J a ( x ~y) f(y) d y- 

If / G I° + (L p (a,b)) with < a < 1 and p > 1 then the left-sided Riemann-Liouville 
derivative of / of order a defined by 

for almost all cc G (a,b). 
For G (0, 1), the integral transform 



(ifH/)(t) = f K H (t,s)f(s)ds 
Jo 
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H+- 

is a isomorphism from L 2 ([0, 1]) onto I Q+ 2 (L 2 ([0, 1])) and its inverse operator Kjj 
is given by 



K H l f = t H ^D* + H l 2- H f> for H > 1/2, 
K H l f = t l 2- H Dl; H t H - L 2D™f for H< 1/2. 



(2.4) 
(2.5) 



Moreover, for H < ^, if / is an absolutely continuous function then K^f can be 
represented of the form ( see [9] ) 



(2.6) 



3 The maximum likelihood estimator and Cramer-Rao 
type bound 

We consider a function 9 = (0 1 , . . . , 6 d ^j belonging to M. An estimator £ = . . . , 
of = (0 1 , . . . , 6 d ) is called unbiased if, for every t G [0, T] 

Ee(£) = 6l i = l,...,d 

and it is called adapted if, for each i = 1, . . . , d, £ l is adapted to (^7) tg j r j. 
Since for any i = 1, d, the function 0* is deterministic and 



Jo 



(K H L (6 l )(s)yds < oo, 



then Girsanov theorem yields that there exists a probability measure Pg absolutely 
continuous with respect to P under which the process B H := (b^; t € [0, T]j defined 
by 



te[0,T) 



(3.7) 



is a d-dimensional fBm with Hurst parameter H and mean zero. Moreover the Gir- 
sanov density of Pg with respect to P is given by: 



dPe 
dP 



exp 



W f 2 K^MdWi 



{K^^mfds 
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and 



B? = f K H (t,s)dW s 
Jo 



'0 

where W is a d-dimensional Brownian motion under the probability Pg and 

Wt = Wt- I KgWWda, i = l,...,d; t e [0,T\. 
Jo 

The equation (|3.7p implies that B H is an unbiased and adapted estimator of 9 under 
probability Pg. In addition, we obtain the Cramer- Rao type bound: 

r-T _ rT rp2H+l 

R(H,9) :=K(9,B H ) = E e \\B t H fdt = d t 2H dt = — — -d. 

Jo Jo 2H + 1 

The first main result of this section is given by the following proposition which 
asserts that 9 = B H is an efficient estimator of 9. 



Theorem 1 Assume that H < \. Ift is an unbiased and adapted estimator of 9, 
then 

Ee / Ut-9 t \\ 2 dt>R(H,9). (3.8) 
Jo 

Proof: Since £ is unbiased, then for every ip G M we have 

E v (%) = E v (<4), j = l,...,d. 

Let (p = 9 + eip with ip G M and e G R. Then for every t G [0, T] and j G {1, . . . , d}, 
we have 

= E e+£ ^9l+e4) 
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This implies that for every j = 1, . . . , d 



d 



de/ £ =o 
' d 



E 



de / e =o 



exp 



W f K-\V + e^){s)dWl 
i=i ^° 



- j {K^i? + e^){s)Y)ds 



x(e t -0i) 



vi=l 



Applying Cauchy-Schwarz inequality in L 2 (f2, dPg), we obtain that for every t 6 [0, T] 

j= i j= i • \\- J o 

j= i L Jo J 

We take for each j = 1, . . . ,d, ip 3 t = t 2H . Since t — ► t 2H is absolutely continuous 



function, then by (|2.6p . a simple calculation shows that 



K H \t 



2H(5{ 1 - H,H + 2 ) .H-i/2 
r(± - fl) 

2i?(r(i + H))t H - 1 / 2 . 



It is known that 



< < 1 for every z € [1, 2]. 



(3.9) 
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Combining the facts that zF(z) = T(z + 1), z > 0, 2H < (H + \) 2 and (|3"U]) . we 
obtain 

^ = ||^|| 2 < (r(^+ir)) 2 £; (||et-0 t || 2 )<^(||e t -^ll 2 )- 

Hence, by an integration with respect to (it, we get 

rp2H+l r T 



rp2M+L rl 

R(H,9) = ^— - < E e J Ut - t fdt. 



o 

Therefore (13.8j) is satisfied. 

Corollary 1 The process 9 = B H is a maximum likelihood estimator of 9 . 

Proof: We have for every ip S M 
d 



d 

T , ex P 

de / e =o 



£ f K-^ + s^dWi- 1 - f + e^)(s)) 2 )ds 
l=1 Jo 2 Jo 



0. 



Hence 



£Qf K H \n S )dwi- jT i^^x^^xs)^) =o. 

Which implies that for every i = 1, ...,d 

E H* K^mWdWj - ^Kx\r)(s)Kx l {9%s)ds^ = 0. 

Then, for each i = l, ...,d 

Wt= f K H \9%s)ds, te[0,T\. 
Jo 

Therefore by (|2.3p . we obtain that B = 9. 

4 Superefficient James-Stein type estimators 

The aim of this section is to construct superefficient estimators of 9 which domi- 
nate, under the usual quadratic risk, the natural MLE estimator B . The class of 
estimators considered here are of the form 

S(B H ) t = Bi I + g(B t H ,t), te[0,T] (4.10) 
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where g : M. d+1 — ► R rf is a function. The problem turns to find sufficient conditions 
on g such that 1Z (9, S(B H )) < oo and the risk difference is negative, i.e. 

AK(6) = U (9, 5{B H )) - n (9, B H ) < 0. 
In the sequel we assume that the function g satisfies the following assumption: 



(A) 



J, 1 \\g(B t H ,t)\\ 2 d dt 



< oo, 



the partial derivatives dig 1 := i = 1, . . . , n of g exist. 



Then K (9, 5(B H )) < oo. Moreover 

r r T 1 

AK(0) = E e / \\B? + g(B?,t)-d t \\ 2 d -\\B?-e t \\ldt 
Jo 

\\g{B»Ml + ^(9\B»MB^-ei)) dt 



In addition, 



Ee 



J ° i=l v ' 



d r T 



i=i Jo 



g\x l ,...,x*,t)(x*-9l) 

£^l(^-9j)2 \ 

x e 2t2« dx 1 . . . dx d 



V f T (2^ H )-i ( [ t™d t g\xi 
~{ Jo \Jr* 



, . . . , x , i) 



x e *? H dx 1 . . . dx d (it 



E / 3 {t 2H E e d i g i (B?,t))dt = E e [E I* t 2H d l9 \B^t)dt 

i=i - 70 L<=i - 70 



Consequently, the risk difference equals 



o 



|| 5 ( J Bf,t)|| 2 + 2t 2i? E^(5f,t) U 



i=i 



(4.11) 
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We can now state the following theorem. 

Theorem 2 Let g : R rf+1 — ► M. d be a function satisfying (A). A sufficient conditions 
for the estimator (B^ + g(B^ !^)) te [oT] ^° dominate B H under the usual quadratic 
risk is 

E e £ (\\g(Bf^ 2 + 2f H J^d i g i (Bf,t)\dt < 0. 
As an application, take g of the form 



g(x, t) = at 



2H 



r \\x\ 



(4.12) 



where a is a constant strictly positive and r : R + — > M + is bounded derivable function. 
The next lemma give a sufficient condition for g in (|4.12)) to satisfies the assumption 
(A). 

Lemma 1 If d> 3 and H < h then 



E 



l^ll 2 



dt 



< oo. 



(4.13) 



Proof: Firstly the integral given by (|4. 13[) is well defined, because 

(dt x P) ((t,w);Bl I (w) = 0) = 

where (dt x P) is the product measure. 

Using the change of variable and d > 3 we see that 



E 



1 



\B 



ff 112 



dt 



dt 

~2H 



t 



Hall 2 



2iT\\y\ 



-dy<C I -pijdt, 



where C is a constant depending only on d. Furthermore, since H < \ then (14TT3D 
holds. 



Theorem 3 Assume that d> 3. If the function r, the constant a and the parameter 
H satisfy: 

i) < r(.) < 1 

ii) r(-) is differentiable and increasing 
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Hi) < a < 2{d - 2) and H < 1/2, 
then the estimator 

5{B H ) = Bf - at™ r -10p- B», t e [0,T\. 

dominates B . 



Proof: It suffices to prove that A1Z(8) < 0. From (|4.1ip and the hypothesis i) and 
ii) we can write 



aEa 



^ f '^(|gg!)_ 2(d _ 2) Kii^ii 2 ) 



ir'(\\B?f))dt 



< a [a - 2{d - 2)] E e 



a / t 



r(\\B( 



H\\2\l 



\B?\\ 2 



Combining this fact with the assumption Hi) yields that the risk difference is negative. 

Which proves the desired result. 

For r = 1, we obtain a James-Stein type estimator: 

Corollary 2 Let d > 3, < H < ^ and < a < 2(d — 2). Then the estimator 



at 



2H 



\B?\\ 2 



H 



te[0,T] 



dominates B . 
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